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tj^ ■ Abstract 

To distinguish between purely fractionally integrated (FI) processes, we propose in this article 
a new and appropriate fractional Dickey-Fuller (F — DF) test. This new test extends the familiar 
Dickey-Fuller (1979) type tests for unit root (1(1) against 1(0)) by embedding the case d — and 
d = 1 in continuum of memory properties. The new F-DF test is easy to implement and is based 
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on two time domains properties of FI(d) processes. First, if a time series follows FI(d) process, 
than the (— l + d)th differenced series follows an 1(1) process. Second, the purely FI(d) processes 



are non explosive for d E R. These two properties will allow us to draw a bridge between the 



process 1(0) and 1(1) by testing the general hypotheses test H : d > d against H : d < d , with 
]-§,§[ and d e [0,1]. 

Keywords: First autoregression, fractional integration, explosive processes, unit root, fractional unit 
root. 



Introduction 

Recently, extensively worked area in time series analysis has concerned to the design of appro- 
priate test statistics to distinguish between 1(0), 1(d) d £ ]0, 1[ and 1(1) behavior. Some recent 
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contributions on this topic include Dolado, Gonzalo and Mayoral (2002), Nilsen and Johansen 
(2010), Lobato and Velasco (2007). The fractional unit root distributions was first considered by 
Sowell (1990), who analyzed the behavior of the usual Dickey-Fuller type regression, when the 
errors are fractional. Specifically, Sowell (1990) considered the regression 

Vt = Vt-i + e t , for t = 1, 2, • • • , n, 

where y = and e t is a stationary fractionally integrated process that is, 1(5) process with 
—0.5 < 5 < 0.5. He showed that the ordinary least squares (OLS) estimator for (coefficient on 
yt-i, hereafter </> n ) has a nonzero density over the whole real line for the special case of 5 = 0, i.e., 
a unit root. For other values of 5, <p n has different normalized constant depending on the magnitude 
of 5, i.e. n mm ^' 1+2 ^((f) n — 1) = O p (l) and his asymptotic is function of two distributions which 
both depend on fractional Brownian motion. Furthermore, the standard t statistic of <j) n (hereafter 
t^ n ) converges to the well defined density when 5 = 0, for other values of 5, the asymptotic 
distribution of t-z n diverges to infinity. Sowell (1990) concluded that if his distributional theory 
is used to test the presence of unit roots in fractional ARMA models the implementation would 
requires tabulations of the percentiles of fractional Brownian motion conditionally on 5, for the 
statistic 72 mm [ 1 . 1 + 25 l (0 n — i) 5 and thus might suffer from misspecification. He concluded, also,he 
statistic t2 are not useful. 

Diebold and Rudebush (1991) examined the properties of Dickey-Fuller test under fractionally 
integrated alternatives and showed by Monte Carlo simulations that this test has quite low power 
and can lead to the incorrect conclusion that a time series has a unit root also when this is not true. 
They pointed out that a more appropriate testing procedure is needed to draw conclusions about 
the presence of the unit root. The lack of power of Dickey-Fuller test to distinguish between the 
1(1) null hypothesis and the fractional integrated of order d (FI(d)) alternative, has motivated the 
development of new testing procedures that take this type of alternative explicitly into consideration. 
One approach for testing against fractional alternatives belongs to the Lagrange-multiplier (LM) 



framework studied in Robinson (1991, 1994), Agiakloglou and Newbold (1994), Tanaka (1999), 
Breitung and Hassler (2002) and Nielsen (2004). The aim of this approach is adapting time domain 
procedures and embedding the models of interest in general FI(d) framework, instead of the 
autoregressive alternatives typically considered in the literature. Robinson (1994) has developed a 
test for unit roots, that unlike the familiar Dickey-Fuller test, was not embedded in AR structures 
of form: 

(1-^)^ = ^,4 = 1,2,..., (1.1) 

where L is the lag operator (i.e. Lx t = x t _i) and u t is stationary 1(0) process, and with, the unit 
root null corresponds to: 

H : cf> = 1, 
but based on fractional alternatives of form: 

A d x t = u t , 4 = 1,2,..., (1.2) 

where d can be any real number, A d = (1 — L) d is the fractional difference operator defined by its 
Maclaurin series (by its binomial expansion if d is an integer): 

(-] _ T\d _ V^oo r(-d+j) 

K 1 ^) 2^=or(-d)r(j+i)> 

where 

| f n +0 ° s'^e-'ds, if*>0, 

r(,) = 

I oo if z = 0, 

if z < 0, T (z) is defined in terms of the above expressions and the recurrence formula zT (z) = 
r (z + 1) and u t is stationary 1(0) process, and where the unit root null corresponds to: 

H :d=l. 

Robinson (1994) has also proposed a Lagrange Multiplier (LM) test where the null hypothesis is: 

Hq : d = do 
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in model given by (1.2) for any given real value d . 

The major consequence of this change, in the set of alternative models considered, is that the 
the asymptotic distributions are standard. Nonetheless, the advantage of having a standard limit 
distribution, Tanaka (1999), showed with simulation experiments that the LM tests have the serious 
size distortion (see Tanaka 1999 fore more detail). Another critic addressed for the LM tests is 
that, by working under the null hypothesis, it does yield any direct information about the correct 
long-memory parameter d, when the null is rejected (Dolado, Gonzalo and Mayoral 2002, Bertrand 
candelon, Gil Alana 2003). In order to overcome the drawback, researches are directed towards a 
global procedure test, which embedded the AR structure form (1.1) and the fractional form (1.2) 
in the same model. The aim of this approach is to extends the well known Dickey and Fuller 
approach, originally designed for the H : d = 1 against H : d = 0, to the more general setup 
H : d = d against H : d < do, with < d < 1. 

The fractional Dickey-Fuller (FD-F) unit root test was first considered by Dolado, Gonzalo and 
Mayoral (2002), (hereafter DGM). The (DGM) test for the null hypothesis d = d against a simple 
alternative d = d\, with d\ < d , is based on the OLS estimation of the following model 

A do x t = pA dl x t - 1 + u t , (1.3) 

where u t is 1(0) stationary process. When d = 1, the model (1.3) becomes 

Ax t = pA dl x t -! + u t , 

where d\ < 1. When di is not taken to be known a priory, a pre-estimation of it is needed to 
implement the test. Lobato and Velasco (2006) show that the DGM test is inefficient and the 
regression model (1.3) is misspecified because it does not include the data generating process 
defined by (1.1) and (1.2) like a particular case. The regression equation, proposed by DGM, is 
used only to suggest a test statistic, and we can not considered it as a generalization of familiar 
Dickey-Fuller test for unit root. 



To design the appropriate test statistics to distinguish between 1(d) process, we propose in 
this article more flexible and more appropriate regression model for fractional unit root test. By 
using the non explosive feature of ARFIMA(0,d,0) processes, our testing procedure generalize the 
familiar unit root test in the most adequate way. Indeed, there is a significant difference between 
the autoregressive (AR) models of form (1.1) and the fractional alternatives of form (1.2). As noted 
by Gil-Alana and Robinson (1997) and Gil-Alana (2004), "fractional departures from (1.1) and 
(1.2) have very different long run implication. In (1.1) for \<J)\ > 1, x t is explosive, for \<f)\ < 1 
x t is covariance stationary, and for = 1 it is nonstationary but not explosive. In (1.2), x t is 
nonstationary but non explosive for all d > 0.5. As d decreases beyond 0.5 and through 1, x t can 
be viewed as becoming "more nonstationary" (in the sense, for example, that the variance of partial 
sums increases in magnitude) but it does so gradually, unlike in case of (1.1) around <\> = 1." 

The non explosive feature of ARFIMA(0, d, 0), often quoted in the literature, has neither been 
studied theoretically nor been used in practice, for the needs in statistical inference on the fractional 
unit root test. Intuitively, the non explosive feature of the process ARFIMA(0, d, 0) for all d > 0.5 
means that if a first order autoregression model is fitted to a sample of size n generated according 
(1.2) then the OLS estimator of the first order autoregression parameter will not exceed 1. Figure 
1 below illustrates this fact in an obvious way. 
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This figure was made as follows: For a fixed sample {«!_„, • • • ,u , ■ ■ ■ ,u n } generated from 
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i.i.d.N(0, 1), with n = 1000, samples of ARFIMA(0, d, 0) processes were generated for d varying 
between 0.5 and 2, with step of 0.01. For each sample {x t , t — 1, • • • , n} a first order autoregression 
model 

x t = $ n x t -i + % (1.5) 

is fitted and estimate of (hereafter 0„) are calculated. By plotting the parameter <f> n against the 
fractional differencing parameter d, one obtains the figure I. This figure indicates that if the value 
of the parameter d is lower than I, then the value of the parameter <p n is lower than I. While, when 
the parameter d is greater or equal to one, the value of the parameter <p n converges to I. Based 
on the relation between <fi n and d highlighted by figure 1, the non explosive feature of the process 
ARFIMA(0, d, 0) will allow us to build a unified theoretical framework for the unit root test and 
fractional unit root test. More precisely, to distinguish between FI(d), d G ] — \, § [ processes, we 
present a new and appropriate Fractional Dickey-Fuller test in time domain that extends the familiar 
Dickey-Fuller (1979) type tests for unit root (1(1) against 1(0)), by embedding the d — and d — 1 
in continuum of memory properties. The proposed statistical test is based on the following idea: 
to test if a given process y t is integrated of order d (y t ~» 1(d)), it is enough to test if the process 
x t — (1 — L)~ 1+d y t is integrated of order 1. This simple idea associated with the non explosive 
feature of the ARFIMA(0, d, 0) processes will allow us to draw up a bridge between the process 
1(0) and 1(1), by using a new fractional autoregression model (see section 3), which includes 
process (1.1) and (1.2) like particular case and which makes it possible to test the following 
general testing problem 

H : d > d against H± : d < d , (1.6) 

with d e ] -|, | [ and d e [0, 1]. 

The rest of the paper is organized as follows. The next section contains an introduction to some 
standard results and concepts of both fractional integrated series and the fractional space D [0, 1], 
that we need for this study. In section 3, to explain the relation between n and d, highlighted by 
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figure 1, we will study the behavior of the asymptotic distribution of n when d G [|, oo[. In the 



same section, we discuss the consequences of the nonexplosiveness of ARFIMA(0, d, 0) process 
to deal with the unit root test and fractional unit root test. In section 4, by using a new auxiliary 
fractional autoregression model, we provide the main results on asymptotic null and alternative 
distributions for the testing problem (1.6). In section 5, through Monte-Carlo study, we show that 
the proposed test fare very well both in terms of power and size when we use the t— statistic. In 
section 6, we give some concluding remarks. 



In this section, the expression D [0, 1] denotes the space of functions on [0, 1] in which all 
elements are right continuous and have left-hand limits, endowed with the Skorohod topology (see 
Billingsley, 1968, p. 111). Weak convergence of probability measures on D [0, 1] and convergence 



The main technical tool, that we need for this study is in Sowell (1990) and Liu (1998). To 
begin with, we briefly outline some results, which will be used extensively, concerning the limit 



and u t: t = 0, ±1, • • • , are i.i.d. random variables with E (u t ) = and E\u t \ a < oo for a > 
max {4, -85/ (1 + 25)}. Since 



I. Preliminaries 




behavior of sample moments of long memory processes. Let y t = (1 — L) 5 u t , with 5 G ]— 0.5, 0.5] 




for S = 0.5 



for - 0.5 < 5 < 0.5 



(2.1) 



from the invariance principle of Davydov (1970) it follows, as n — > oo 




(2.2) 



where 




(2.3) 



and B 5)0 (r) is the type I fractional Brownian motion on D [0, 1] defined as follows, 

ajjy j ^ Ut - sf - (s) 5 } dB(s) + / r (t - s) s dB(s) if 8 e } -I \ [ , 



B 5 ,o(r) 

rB(l) 

where B(s) is standard Brownian motion and 



if 5 = 0.5 



1/2 



For the time series y t — (1 — L)~ s ~ m u t , with m > 1 and 5 e ]— |, §] by combining the 
continuous mapping theorem, the results from Sowell (1990) and Theorem 2.2 of Liu (1998) we 
have the following useful results 



K(5,n)n-V*+ s +™ yinr] ^ B<5 > m(r) ' 



B 5j0 (ri)dr 2 dr 3 • • • dr m , 



K 2 (5,n)n 2 ( m+<5 ) ^ 

v ' 7 = 1 



[B 5>m (s)] 2 rfs, 



(2.5) 
(2.6) 

(2.7) 



where K 2 (5,n) is defined by (2.3) and 



B 5 , m (r) = 



B.nfr) 



><5,0l 

rr rr m -i 



JO 



if m — 1 

J r2 B ( 5 i0 (ri)dridr 2 dr3 • • • dr m _i if m > 2. 



(2.8) 



II. Consequences of nonexplosiveness of the ARFIMA(0, d, 0) processes 
Consider a purely fractionally integrated process {y t } defined by 



1 - L) m+S y t = u t , for t = 1, 2, • • • n, 



(3.1) 



with y = 0, and where m > and 5 e ] — \, |], L is a lag operator and u t is defined as in the 
section 2. We denote m + 5 = d. 

Theorem 1. Let {y t } satisfy (3.1). If the first autoregression model (1.4) is fitted to a sample of 
size n then 



(logn) (J n - l) 



-2 



Jo [ B i o(0 



2 ' 



(3.2) 



when 5 = 0.5 and m — 0. 



n l+25 



_ ri _i_ ai r ( 1+<5 ) 

V 7 Jo B 8, (^ dr 



when —0.5 < 5 < and m—1. 

x ±{B 2 (1)-1} 



i/ 5 = and m—1. 



fc- 1 )^ r4? m r (3 - 5) 



i/ (m = 1 and < 5 < 0.5) or (-0.5 < 5 < 0.5 and m > 2). Here £*, m (r) w de/med £y (2.8). 
Proof. See Appendix 

Theorem 1 indicates that the least squares estimate is super consistent, for d > 1. The rate 
of convergence depends on the order of integration. It is well known that if x t is 1(0), the OLS 
estimate converges at the rate n 1 / 2 , and if x t is 1(1) i.e. 5 = and m—1, convergence is at the rate 
n. This may suggest that the rate of convergence increases with the order of integration. However, 
this is not the case. If a series is I(m + 5) for m > 1 and < 5 < 1/2, <p n converges at the rate n, 
the same for all m and 5. If a series is 1(1 + 5) for — 1/2 < 5 < 0, 4> n converges at the rate n 1+25 . 
So for —1/2 < 5 < —1/4 the rate of convergence is slower than n 1 / 2 . The figure 2 illustrates this 
fact clearly. Furtheremore, when m — 1, the rate at which 7l mm [ 1 > 1 + 2<5 ] (^p n _ converges to its 
limiting distribution is slow for nonpositive values of 5. This implies that n converges very slowly 
towards 1 for d < 1. Moreover, for < 1 the limiting distribution of (4> n — lj has nonpositive 
support and then 

limP U n - 1< O) = 1. 

When d > 1, 7T, mm [ 1 > 1 + 2 ' 5 ] — converges to its limiting distribution at a faster rate than they 
do in standard first autoregressions with stationary variables. This implies that n has high speed 
convergence towards 1 for d > 1. Moreover, for d > 1 the limiting distribution of — lj has 
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nonnegative support and then 



lim P 

n— >oo 



min( I . I -+2&) 



- 1 < 



0. 
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The relation between the first autoregressive parameter 0„ and the fractional parameter d, high- 
lighted by the results of Theorem 1 and illustrated by the figure 1, in the context of ARFIMA(0, d, 0) 
process, suggests that when we deal with unit root test, we have 4> = 1 or <\> < 1, according to 
wether d > 1 or d < 1. In other words, the testing problem 



Hq : (p = 1 against Hi : <p < 1 



(3.6) 



is equivalent to the testing problem 



H : d > 1 against Hi : d < 1 



(3.7) 



It is obvious that the classical testing hypothesis of the Dickey-Fuller test is (3.6). If one would 
test the hypotheses (3.6), then in the familiar Dickey-Fuller test we have 

Hq : (p = 1 implies d — 1, 
Hi : (j) < 1 implies d = 0, 

whereas in the context of ARFIMA(0, d, 0) processes we have 

Hq : = 1 implies d > 1, 
Hi : (p < 1 implies d < 1. 
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Now, we can expose our simple idea associated with the non explosive feature of the ARFIMA(0, d, 0) 
process which will allow us to draw up a bridge between the process 1(0) and 1(1), by using a 
new fractional autoregression model. The new fractional model includes process (1.1) and (1.2) 
like a particular case. For the testing problem (1.6), we propose to test the null hypothesis by 
means of the t-statistic of the coefficient of (1 — L)~ 1+do y t ^i in the ordinary least squares (OLS) 
autoregression 

A- 1+do y t = 0A- 1+ S*-i + € U t = 1, 2, • • • , n, (3.8) 

or equivalently 

A d °y t = P A- l+do y t ^ + e t , t = 1, 2, • • • , n, (3.9) 

where p = — 1 and e t the residuals. If d — d and e t = u t (recall that d is the true value of 
integration parameter and d is the value specified under the null), in the fractional autoregressive 
model (3.8) and (3.9), according to wether 

d — and 0=1 i.e. (p = 0) 
d e ]0, 1[ and = 1 i.e. (p = 0) 

d — 1 and 0=1 i.e. (p = 0) 

we have respectively 

y t = u t i.e. y t ~* 1(0) 

(1 - L) d °y t = u t , d e ]0, 1[ i.e. y t - FI(d ) 

(1 - L)y t = u t i.e. y t -w 1(1) 
In order to grasp the intuition behind the fractional autoregressive model (3.9), suppose that 

y t FI(d ) and let us consider the relation between A d °y t and A -1+do y t _ 1 . Note that, it is easy 

to check, that A da y t = (1 - L) [A- 1+da y t ] and A- 1+d "y t ~* 1(1). Putting A- 1+da y t = x t we can 

rewrite (3.9) as follows 

(1 - L)x t = px t -i + e t , t = 1, 2, • • • , n. (3.10) 
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The regression model (3.10) is the simple Dickey-Fuller framework to deal with the testing problem 
(3.6). Furthermore, suppose that y t -w FI(d) then x t ~* FI(1 + d — d ) and we have the following 
(because of nonexplosiveness of ARFIMA(0, d, 0) processes) 

= 1 i.e. p = 0, 

(f) < 1 i.e. p < 0, 

according to wether 

1 + d — d > 1 i.e. d > d , 

1 + d — d < 1 i.e. d < d . 
In other words, by using the regression model (3.9), the testing problem (1.6) is equivalent to 

the testing problem (3.6). 

The simple fractional autoregression model (3.9) can be easily implemented for practical settings 

and is flexible enough to account for broad family of long memory specification of the fractional 

parameter d. The OLS estimator of p (hereafter p n and its t -ratio (hereafter tp n ) for the regression 

model (3.9), are given by the usual squares expressions 

~ _ EILi [A*yJ [A-^yt-i] 

P» r -i 1/2 ' 

{^ELi [A- 1+ *y t -i] 2 } 
where the variance of the residuals, s 2 n is given by 

2 E?=i l A "°^ - pA' 1+d "yt-i] 
S ™" n 

We call the testing procedure based both on the regression model (3.9) and the hypothesis test 
(1.6) the "Fractional Dickey and Fuller" (F - DF) test. 

III. Asymptotic null and alternative distribution and consistency of F-DF test 

A. Asymptotic null and alternative distribution 

The asymptotic distribution of the statistic's p n (appropriately standardized) and tp n , under the 
hypothesis H : d > d and if i : d < d , is given in the following theorems. 
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Theorem 2. Let y t satisfy (3.1). For 5 = d — d e ] — |, | [, /f^e regression model (3.9) is fitted 
to a sample of size n then, as n — > oo 

n 1+2 % =► -^4rM » if - 0.5 < 5 < 0, (4.1) 
Jo B 2 (r)dr 

i{B 2 (l) - 1} 

^P»^ f i ' \ , » if5 = > (4-2) 
Jo B 2 (r)rfr 

n P» h L ~oV7 J » if 0.5 > 5 > 0, (4.3) 



1 {B 2 (l)| 
/oB2(r)dr' 



where B$(r) and B(r) are respectively the type I fractional Brownian Motion and the standard 
Brownian Motion. 

Theorem 2 indicates that the rate of convergence depends on the difference between the unknown 
true value of the parameter d and the value specified under the null d i.e. 5 = d — d . The rate of 
convergence increases with 5 for — | < 5 < until the rate n, and converges at the same rate n for 
all 5 > (See figure 2 above). The asymptotic distributions of p n depends also on the 5 = d — d . 
If S = i.e. (d = d ), then the asymptotic distribution of p n is reduced to that deduced by Dickey 
and Fuller (1979) for the particular case H : d = 1. In the case 5^0, then the distribution of p n 
are reduced to those deduced by Sowell (1990) for particular case H : d — 1 and d — 1 = 5 with 

Theorem 3. Let y t satisfy (3.1). For 5 = d — d e ] — |, | [, regression model (3.9) is fitted 
to a sample of size n then, as n — >■ oo 

tp n A -oo, if - 0.5 < 5 < 0, (4-5) 

^iS^' if5 = ' (46) 

[Jo B 2 (r)drj 

A +oo, if 0.5 > 5 > 0, ( 4 - 7 ) 

where B(r) is the standard Brownian Motion. 

Theorem 3 indicates that the asymptotic distribution of tp n depends also on 5 = d — d . If 5 = 
the asymptotic distribution of tp n is reduced to that deduced by Dickey and Fuller (1979) for the 
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particular case H : d = 1. In the case 5 > and 5 < the distribution of tp n diverge respectively 
to (+oo) and (— oo). Similar results, for the particular case H : d = 1 with d — 1 = 5^0 are 
derived by Sowell (1990). 

Corollary. Let {y t } be generated by (3.1). For 8 — d — d e ] — \, \ [. If regression model (3.9) 
is fitted to a sample of size n then, as n — >■ oo, we have for the t-statistic 



and 



{tpn) =>■ -<^6 , r i TT72 if 5 < 0, 



where = 



r(i-2<5) 



1/2 



and where Bg(r) and B(r) are respectively the type I fractional 



_(i+25)r(i+5)r(i-(S) 
Brownian Motion and the standard Brownian Motion. 

The limiting distribution of rH 5 ! (t^ n ) has nonnegative support if 5 > and a nonpositive support 

if 5 < 0. 

The standardized least squares estimate p n and the corresponding ^-statistic, tp n will be noted 
hereafter respectively by r p and r t . By using these notations, we explain now, how one can exploit 
the results of theorem (3) to statistical inference in time series. Our interest is to show that the 
F-DF tests, based on r t statistics are consistent whereas the F-DF tests based on r p statistics are 
not consistent. 



B. Consistency of F-DF test based on r t statistic 

Consider the problem of hypothesis test (1.6) in sample of size n. It is convenient to introduce 
the nonrandomized test defined by a function <fr n on the sample space of the observations r t , with 
critical region C. The test for a region C is its indicator function 

l if neC 

®n (n) = < (4.8) 

o if r t ic 
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Let a and (5 respectively the type / error and the type 77 error of the test <fr n . Since H and Hi 
are composite, we have 

a = Sup d > d() n (d) 
(3 = Inf d<do 'Ki(d) 

where ir (d) = P Ho (reject H ) and n 1 (d) = P Hl (accept H ). Practically, this entails computing 
a statistic r t from a sample, whose distribution P Ho when the null hypothesis H is true can be 
tabulated, and used to fix the probability of rejection when H is true (a Type I error) not exceed a 
chosen value a. Since, the alternative hypothesis is H 1 : d < d , it is natural to consider one sided 
critical regions of the form 

C = {r t <c n (a)}, (4.9) 

where a is the level of the test. With these settings, the power function of the test <fr n , denoted by 

7r#„(d) is : 

7r*„(d) = 1 - Infa^n^d). 

Theorem 4. A sequence of tests {$> n } defined by (4.8), with critical region (4.9) each of given 
size a is consistent i.e. 

lim 7r$> n (d) — > 1 for d < d . 

n— s-oo 

Proof. If 5 < 0, from theorem 3 and corollary 1 it follows 

lim P Hl (r t < c n (a)) = lim P Hl (n 5 T t < n s c n (a)) 

n — >co n — >co 

= iim p Hl (~cs < o] 

= 1, because n s r t has nonpositive support. 

As n becomes infinite, ^^ n {d) tends to 1, if 5 < 0. Moreover, for each n, n® n {d) is an increasing 
function, with larger absolute value on 5 < than 5 = 0. Hence, for testing d > d against d < d , 
the region r t < c n (a) is unbiased (for any n) and defines a consistent family of tests. 
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We have also 

lim P Ho (r t > c n (a)) = 1 - a, if 5 = 0, 

n — >oo 

and 

lim P Ho (r t > c n (a)) = 1, if 5 > 0. 

n — >oo 

For the case where 5 = 0, the exact critical points c n (a) are given by Fuller (1976, pp. 373 
and 375). These critical points are also asymptotically valid (with in the case where the u t 's are 
not normal). This implies that the proposed test can be understood and implemented exactly as the 
Dickey and Fuller test for unit root by using the usual table statistics. 

In practice, the statistician or econometrician, would wish to know the answer to the following 
question : Is the series, under study, stationary or not ? To answer this question, we propose the 
following downward testing procedures where the highest integration level is tested first. 

1) Reject H : d > 1 (against H 1 : d < 1) and go to step 2 if r t < c n (a); otherwise, conclude 
that H : d > 1 is true. 

2) Reject H : d > 0.5 (against Hi : d < 0.5) if r t < c n (a); otherwise, conclude that Hi : d > 
0.5 is true and 0.5 < d < 1. 

IV. Power of F-DF test in finite sample 

After the theoretical analysis of the F-DF test based on r t statistic, we now conduct a Monte- 
Carlo experiment to examine the finite sample performance of F-DF test based on r t and the F-DF 
test based on r p . 

We consider the data generating processes, 

(1-L) d y t = u t . (5.1) 
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The regression model will be used for estimation and inference to deal with the hypothesis test 

Ho : d > d against H : d < do, (I) 

is 

(1 - L)*>y t = p n (l - L)~ 1+do y t _i + e t . (5.2) 

Before giving the main results of these experiments, we examine initially, the two following simple 
hypothesis tests 

Ho : d = d against Hi : d < d (II) 

and 

Ho : d = do against Hi : d > d (HI) 

In the experiments reported in this section, a general procedure for generating a stationary frac- 
tionally integrated series of length n is to apply, for t — 1, • • • ,n and some fixed m, the formula 

Xt = g r ( «oro- + i)"-* (5 - 3) 

where {ui_ m , ■ ■ ■ ,u n } is a random sequence of i.i.d.N(0, 1). This simulation strategic for type I 
fractional Brownian processes has suggested by Davidson and Hashimzade (2008). They argued 
that, by choosing m large we should be able to approximate the type I processes to any desired 
degree of accuracy. For the non-stationary parameter confugarations d = 0.6, 0.7, 0.8, 0.9, 1, 1.1, 
1.2, 1.3 and 1.4, we use the following: 

Xt = Xt ~ 1+ g rWTI) ^' (5 ' 4) 

where 5 G ]— \,\[- For the particular case d = 0.5 we can use (5.3) or (5.4) and the value of 
m = 500. All computations were done in Eviews 4.0. 



18 

A. The size of the hypothesis tests (II) and {III). 

Trough Monte Carlo study we show that the F-DF test, based on the fractional first autoregression 
model (5.2), fare very well in terms of size when we use r p or r t statistics. To investigate the size of 
the hypothesis test (II) and (III), 10000 samples of FI(d) Gaussian processes (5.1) are generated 
and the regression model (5.2) is used for estimate r p and r t . The samples sizes considered are 
n = 50; 100; 150; 200 and 250. Samples of n observations were generated for five values of d : 0; 
0.2; 0.5; 0.8; 1 and for each value d, we specify under the null only one value d equal to d. 

For the case where u t 's are normal, the exact critical points for r p and r t , under the null are 
given by Fuller (1976, pp. 373 and 375: Table 8.5.1 for r p and Table 8.5.2 for r t ). For the testing 
problem (II) we consider one-sided critical regions of the form r p < b n (a) and r t < c n (a). For 
the testing problem (///) we consider one sided critical regions of the form t p > b n (l — a) and 
r t > c n (l — a), where a is the level of the test. 

Table 1 reports the acceptation frequencies at, respectively the a = 1%, 5%, 10% significance 
level, for the testing problem (//) when we use the dickey Fuller ^-statistic r t . It is noted that 
the estimated frequencies are very close to the theoretical frequencies (1 — a = 0.99, 0.95, 0.90). 
Similar results, not reported here, are found for the test hypothesis (///) when we use the dickey 
Fuller ^-statistic r t and for the tests hypothesis (//) and (///) when we use the Dickey-Fuller 
'normalized biais' statistic r p . 
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Table 1. 





T\d 





0.2 


0.5 


0.8 


1 




50 


98.96 


98.93 


98.91 


98.90 


99.12 




100 


98.94 


98.91 


98.96 


98.93 


99.12 


a = 0.01 


150 


99.05 


98.98 


98.96 


98.81 


99.04 




200 


99.05 


98.90 


99.11 


99.01 


99.08 




250 


99.01 


98.83 


99.18 


99.11 


98.89 




50 


94.05 


95.22 


95.22 


95.10 


95.13 




100 


94.48 


95.05 


94.97 


94.95 


95.17 


a = 0.05 


150 


95.09 


95.09 


94.79 


94.59 


94.91 




200 


95.04 


95.20 


95.32 


94.83 


94.88 




250 


94.68 


94.90 


95.44 


94.93 


94.87 




50 


90.34 


89.79 


90.00 


90.25 


90.02 




100 


89.73 


99.76 


90.21 


89.97 


90.24 


a = 0.10 


150 


90.17 


90.34 


89.86 


89.72 


89.95 




200 


90.39 


90.23 


89.84 


89.66 


89.99 




250 


89.49 


90.04 


90.49 


90.11 


90.15 



Although these results are very significant we have also, used the samples of 10000 observations 
to estimate the densities (following Sowell (1990)) of t p and r t under 5 = by the kernel estimator 

1 10000 / _ ~ \ 

^) = Toooo5g*( -IT )< 

where 6j are the estimated values of t p under (5 = 0) or r t under (5 = 0) for the 10000 samples 
and ^ (y) is the density defined by (i§) ^(1 — y 2 ) 2 ^j for — 1 < y < 1 and zero elsewhere. The 
value of h was chosen to minimize the integrated mean square error (see Tapia and Thompson 
(1978, pp 67)). The estimated densities are presented in Fig. 3. and Fig. 4. 

For each statistics r p or r t under 5 = (i.e. d = d ) and for a given size of sample n, the 
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estimated densities for different values of d are represented on the same graph. The figures 3 and 

4 shows that by fitting the regression model (5.2) to the sample generated according (5.1), one 
obtains the same distribution that those used by Dickey-Fuller (1979,1981). 

B. The power of the hypothesis tests (II) and (III) 

To study the power of fractional unit root test, we have generated 10000 samples of FI(d) 
Gaussian process according (5.1) for d = 1. To evaluate the power of the test (II) we specify 
different values of d under the null : 1, 1.1, 1.2, 1.3 and 1.4. The figure below represent the 
different densities of r p under (5 < 0) for the case where the true value of d is equal 1 and 

5 G {—0.4; — 0.3; — 0.2; — 0.1; 0}. Let b n (a) the critical point at the significance level a for r p 
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under (6 < 0), we have 



P [r p > b n (a) \8<0]=P'. 



If the test (II) has good performance in terms of power then the probability (3' should decrease 
when 5 decreases (i.e. d increases) and moves away from zero. f3' should be always less than 
(3. The figure 5 below show the opposite i.e. 0' increases when 5 decreases, because, as noted 
by Sowell (1990), the density appears to be converging to unit mass at zero when 5 approaches 
(—0.5). Consequently, we can not use the statistic r p for the hypotheses test (//). 



Now considering the test (///) in the case where d = d + 5 is the true value of fractional 
integration of the process y t with 5 G {0; 0.1; 0.2; 0.3; 0.4}. The figure below represents the 
estimated densities of r p under (5 > 0) for a size of sample n = 250, if the true value of the 
fractional parameter of integration is equal to 1, whereas the values specified under the null are 
d e {1; 0.9; 0.8; 0.7; 0.6}. The two graphs show that the various densities t p under (5 > 0), with 
5 G {0; 0.1; 0.2; 0.3; 0.4} are all to fit together. In other terms, when Hi : d > d is true, the 
distribution of r p under (5 > 0) does not make it possible to have a powerful test. Contrary to the 
preceding test, the reason of disappointing results of the third test is not in the distribution of r p , 
but in the parameter p. Indeed, the parameter p takes the value if the null is true, i.e. d = d and 
it takes also value if the alternative is true i.e. d > d . 



2(1 



1.6 - 



1.2 - 



OS - 
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C. Sz'ze and power of the F-DF test based on hypothesis test (I) and r t 

In this subsection, in Monte Carlo study we show that the proposed hypotheses test (/) based on 
the DGP (5.1) and in the auxiliary regression model (5.2) fare very well both in terms of power 
and size when we use the r t statistic. To investigate the size and power of the hypotheses test (/), 
10000 samples of FI(d) Gaussian processes (5.1) are generated and the regression model (5.2) 
is used to estimate r t . The sample size considered is n = 50 and n = 250. We will use, as true 
values of the fractional parameter of integration of the process y t , three values d : 0; 0.5; 1 and for 
each value, we specify various values for d . If the set of the values of d for a given value of d 
is denotes by Sd(d ) then the sets which will be used for the three values of d are respectively 

S (d ) = {-0.4; -0.3; -0.2; - 0.1; 0; 0.1; 0.2; 0.3; 0.4} , 
S . 5 (do) = {0; 0.1; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.9} , 
Si (do) = {0.5; 0.6; 0.7; 0.8; 0.9; 1; 1.1; 1.2; 1.3; 1.4}. 

The tables 2 and 3 contains the simulations results on the size of the F-DF test for the hypotheses 
test (/). The tables 2 and 3 shows that the testing problem (/) has good performances in terms of 
size, since we have 

P(r t > c n {a) | 5 > 0) > 1- a. 
The table 5 and Table 6 contains the simulations results on the power of the F-DF test for the 
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hypotheses test (J). There are some conclusions to be drawn from it. First, the power of the F-DF 
test increases with the increase of sample size and 5 = d — d . For examples, for a = 5%, d — 1 
and 5 = -0.1, power is 12.36% for n = 50, 20.76% for n = 250and for a = 5%, d = 1 and 
5 = —0.3, power is 48.5% for n = 50, 86.05% for n = 250. Second, as shown in table 5 and 6, 
for n = 250, the power of F-DF test is below 50% for (5 = -0.1) and for (a = 1%, 5 = -0.2). 
Third, for a given n, a and 5, the power for d = 0, d = 0.5 and d = 1 are approximately similar 
because the asymptotic under the alternative does not depend on d but depend on 5 = d — d . 



Table 2. Size of the hypotheses test (/) when we use r t (n = 50) 







a \ S 


0.4 


0.3 


0.2 


0.1 









1% 


100 


99.99 


99.93 


99.73 


98.96 


do 


= 


5% 


99.98 


99.85 


99.39 


98.17 


94.05 






10% 


99.85 


99.48 


98.26 


95.78 


90.34 






1% 


99.98 


99.97 


99.90 


99.74 


98.91 


do - 


= 0.5 


5% 


99.91 


99.78 


99.28 


98.14 


95.22 






10% 


99.78 


99.42 


98.25 


95.63 


90.00 






1% 


100 


99.99 


99.97 


99.79 


99.12 


do 


= 1 


5% 


99.95 


99.92 


99.55 


98.34 


95.13 






10% 


99.83 


99.54 


98.72 


96.03 


90.02 



24 

Table 3. Size of the hypotheses test (I) when we use r t (n = 250) 
a\ 8 0.4 0.3 0.2 0.1 

1% 100 100 99.99 99.92 99.01 
^ = 5% 100 99.99 99.95 99.29 94.68 
10% 100 99.97 99.73 97.87 89.49 
1% 100 99.98 99.8 99.91 99.18 
d = 0.5 5% ioo 100 99.97 99.27 95.44 
10% 100 100 100 97.66 90.49 
1% 100 100 100 99.92 98.89 
d = l 5% 100 100 99.96 99.32 94.87 
10% 99.98 99.97 99.80 98.00 90.15 



Table 5. power of the hypotheses test (I) when we use r t (n = 50) 

a\6 -0.1 -0.2 -0.3 -0.4 

1% 3.61 9.53 22.23 45.82 

d = 5% 12.66 26.68 47.16 73.03 

10% 21.97 39.56 62.14 84.92 

1% 2.87 9.36 22.59 45.16 

d = 0.5 5% H.9 25.17 48.34 72.74 

10% 21.17 39.58 63.23 84.39 

1% 3.2 9.39 22.91 46.23 

d=l 5% 12.36 26.1 48.50 73.47 

10% 21.63 39.24 63.69 84.99 
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Table 6. Power of the hypotheses test (/) when we use r t (n = 250) 







a \ 5 


-0.1 


-0.2 


-0.3 


-0.4 






1% 


7.59 


30.07 


68.00 


95.48 


d 


= 


5% 


19.98 


52.53 


85.90 


99; 16 






10% 


30.62 


64.77 


92.95 


99.87 






1% 


7.3 


30.06 


68.54 


95.60 


d-- 


= 0.5 


5% 


20.35 


52.40 


85.80 


99.39 






i no/ 
10 To 


O 1 A A 

ol.44 


00. 1U 


no txn 

y/.ou 


nn on 

yy.oy 






1% 


7.71 


29.90 


68.36 


95.39 


d 


= 1 


5% 


20.76 


51.77 


86.05 


99.35 






10% 


31.71 


64.29 


92.41 


99.87 



V. Discussion and concluding Remarks 

In this paper, to distinguish between FI(d) processes with d e ] — \, § [, we have proposed a new 
and appropriate testing procedure in time domain that extends the familiar Dickey-Fuller (1979) 
types tests for unit root (1(1) against 1(0)), by embedding the case d = and d — 1 in continuum 
of memory properties. The main idea of our test procedure is the following: in order to test if the 
process y t is fractionally integrated of order d , it suffices to test if the process x t = (1 — L)~ 1+d °y t 
is integrated of order 1. We have referred to the test based on this original idea as the Fractional 
Dickey-Fuller (FD-F) test. The proposed test is based on the OLS estimator (p n ) or its £ -ratio in 
the autoregression model 

A d »y t = pA- 1+do y t ^ + e t , t = 1, 2, • • • , n 

With this regression model associated with the non explosive feature of FI(d) processes, we have 
showed that the testing problem H : d > d against Hi : d < d with d e ] — \, § [ and d e [0, 1], is 
equivalent to H : p = against Hi : p < 0. We have also, showed that the asymptotic distributions 
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for ordinary least squares (OLS) and its t-ratio under the null simple hypothesis H : d = d 
are identical to those derived by Dickey and Fuller (1979,1981) for the simple case (without drift 
and trend). This implies that the proposed test can be understood and implemented exactly as the 
Dickey-Fuller test for unit root by using the usual tables statistics. It worthnoting that the new 
F-DF test proposed in this paper is a generalization of the particular case H : d = 1 against 
H x : d < 1 (or d > 1) studied by Sowell (1990) for d G ] |, § [, to the general case H : d > d 
against H 1 : d < d with d G ] — |, § [ and d G [0, 1]. For the particular case, Sowell (1990) have 
concluded that the asymptotic distributions of p n and tp n can not be used to test the presence of 
unit root in fractional ARMA models, since the implementation of the test require tabulations of 
the percentiles of fractional Brownian motion conditionally on 5 = d — 1 and thus suffer from 
misspecification. These disappointing conclusions originate from an ill defined statistical problem 
and from an inappropriate use of asymptotic distributional theory. In fact, Sowell have focused 
his attention on the parameter p by considering the test H : p = against H. x '■ P < and 
considered that under the null there are three asymptotic distributions conditionally on 5 = d — 1 
(5 < 0, 5 = 0, 5 > 0). In this paper, by using the non explosiveness of ARFIMA(0, d, 0) processes, 
we have showed, for the general case, that under the null H : p = there are only two possible 
asymptotic distributions conditionally on 5 = d — d (5 = 0, 5 > 0) and under the alternative there 
is only one asymptotic distribution (5 < 0). 

The theoretical framework above, is the unified framework for the unit root test and fractional 
unit root test. Furthermore, in order to test if a given process y t is stationary we can perform the 
downward F-DF testing procedures where the highest integration level is tested first (in our case 
d = 1). We can easily extend our results, by using the same framework, to the sequential testing 
procedure advocate by Dickey and Pantula, allowing the analysis to cover the case d G ] — 1.5, 2.5[ 
and d G [0,2]. 

Further research is currently being undertaken towards generalizing the F-DF testing approach 
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along similar directions as the D-F test has been extended in the unit root literature accounting for 
time series which may exhibit a trending behavior and for general ARFIMA case. 



VI. Appendix 

Proof of Theorem 1. If {x t } satisfy (2.1) and if the first autoregression model x t = <fix t -i +oo t 
or equivalently Ax t = (4> — 1) x t _i +cu t is fitted to a sample of size n then, the least squares slope 
estimate have the following expression: 

l^t=\ X t-\ X^t=l X t-l 

When d = m + 5 = 0.5 i.e. m = 0, 5 = 0.5 and Ax t ~ FI (-0.5) it follows from (2.7), that 

«T 2 (M™ _i y>*-i^ / Bi (r) dr. (A2) 
When m = 1 for -1/2 < 5 < 1/2 and Ax t ~ FI (5), it follows from (2.7), that 

n „i 

k- 2 (5, n) n- 25 - 2 x li / B| (r) dr, (A3) 
t=i Jo 

and when (m > 1 for < 5 < 1/2) and Ax* ~ F7 (m - 1 + 5), it follows from (2.7), that 

k~ 2 (5) n- 2S - 2m x li =► f B| m (r) dr. (A4) 
t=i Jo 

As regards the Ylt=i x t-i^ x t term, we can rewrite it as 

n n n 

\ E W - ^-i - (*« - ^-o 2 } = i E ^» - i E ( Aa ^ 2 • 
<=i *=i <=i 

When m = for 5 = 0.5 and Ax t ~ FI (— 0.5), the first term when multiplied by k~ 2 (5,n) = 
^ logn converges in distribution to \ [B(l)] 2 because relation (2.5), whilst the limiting distribution 
of the second term follows from Lemma 2.1 of Ming Liu (1998) result 2, we have 

n- 1 ^ (A^) 2 (A5) 

71 

When m = 1 for — 1/2 < 5 < and Ax t ~ FI (5), the first term when multiplied by k~ 2 (5, n) n~ l ~ 26 
converges in distribution to \ [B^l)] 2 because of the (2.5), whilst the limiting distribution of the 
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second term follows by using the ergodic theorem 

r 2 (i-5) ' 



n J2t=i (A*,)' 



(A6) 



Therefore, when 5 = 0, i.e., when d = 1, k 2 (5) = a 2 ., B a (r) = B(r) for r e [0, 1] and 

n 

re~ 2 (<J,n) n-^^iA^ | [B 2 (l) - l] , 



(A7) 



whereas when d = m + 5 for (m > 2 and — 1/2 < 5 < 1/2) or (m — 1 and < 5 < 1/2) 

n 

re" 2 (5, n) n 1 " 2 -^ £ x t _ 1 Ax t =► ^B 2 m (1) , 
and when d = 1 + 5 for -1/2 < 5 < 



-2 



(5,n)n 1 y^x t _iAa; f 



^r(i - 25) 
2r 2 (i - 5) ' 



and when d = 0.5 



re 2 (5,n)n 1 Xt-iAxt ==>• —2. 



Hence, using (A2) — (A10) and the continuous mapping theorem, we obtain that 



logn (<j> n - lj 
when d = 0.5. 



re 2 (5,n)n 1 log n Y!t=\ x t-i^ x t 
K~ 2 (5, n) n' 1 J2t=i x t-i 



— 9 



f 1 

Jo 



Bi, (r) 



2 ' 

dr 



1+25 



/2 iA 1 d > 



12 ^ u l r(i-5) 



re 2 (5,n)n 1 E"=i ^t-i^t 



when d = 1 + 5 for -1/2 < 5 < 0. 

re" 2 (5, n) n' 1 EJLi x t-i ^ x t 



n 



(0n - l) 



K- 2 (5,n)n- 2 Eti x Li 
if d = 1, i.e., m = 1 and 5 = 0. 

re- 2 (5,n)n 1 - 2S - 2m J^ =1 x t _ 1 Ax t 



\ [B 2 (1) - 1] 
Jo B2 (0 dr 



n 



(?n - l) 



|B| m (l) 



re- 2 (5, n) n^" 2 ™ £"=i x 2 _, " £ B 2 m (r) dr) ' 
if d = m + 5 for (m > 2 and -1/2 <6< 1/2) or (m = 1 and < 5 < 1/2). 



(A8) 



(A9) 



(A10) 



(All) 



(A 12) 



(A13) 



(A 14) 
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Proof of Theorem 2. The proof is omitted because of its similarity with Theorem 1. Indeed, 
if we denote A~ 1+d "y t = x t then we can rewrite the regression model (3.9) as Ax t = px t -\ + oo t 
with x t ~ FI (1 + 5) and -0.5 < 5 < 0.5. 

Proof of Theorem 3. If y t satisfy (3.1) with d G ] -\, § [ and d - rf = 5 with -0.5 < 5 < 0.5 
and if the regression model (3.9) is fitted to sample of size n then, by denoting A~ 1+d °y t = x t ~ 
FI(1 + 5), the i student statistic have the following expression 

, _ E"=i a:t-i Ax t 



S n {52t=l X t-l} 



1/2" 



We can rewrite as follows 



k(S) 1 n 1 5 £)" =1 x t-i Axt 



1/2 



\k {5)' 1 n-^ {x n f - \k (S)' 1 n- 1 - 5 E^i i^tf 



{k(S) 



- 2 n-^E"=i^ 2 -i} 



1/2 



First notice that 

Y n ]_ n 1 n 2 n 

«n = - (Ax t - p n X t _i) 2 = - (Ax t f + -p^ V] x\_ x p n a; t _iAa; t 

n ^ — ' n z — ' n z — ' n A — ' 

t=i t=i t=i t=i 

Hence, by using A3, A6, A7, A9, Theorem 2 and the continuous mapping theorem, it follows that 



, m-25) 

n r 2 (i-5) u ' 

for all 5 E ]-0.5,0.5[. 

Consider now, the numerator, N, of the tp n statistic, which can rewrite as 



(A15) 



N = 1 -k (6) n 5 [k (Sy 2 n- 1 - 25 [x n f] - 1 -k {d)' 1 
By using (2.5) and (2.7) and A. 7 we obtains 



n 



n 



1 e ( Ax *y 



t=i 



N 



-oo if 5 < 

±{(B(1) 2 -1)} if 5 = 
+oo if 5 > 0. 



(A12) 
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Hence, using A. 3, A. 11 and A. 12 and the continuous mapping theorem, we obtain that 

t ?n 4 -oo if - 0.5 < 5 < 0, 

§{B*(1)-1} 
[Jo B 2 (r)drj 

A +oo if 0.5 > 5 > • 
Proof: Proof of Corollary 1. Direct consequence of Theorem 3. ■ 
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